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Abstract 

We report the first direct observation of the influence of high magnetic fields on the Jahn-Teller 
(JT) transition in TbV04. Contrary to spectroscopic and magnetic methods, X-ray diffraction 
directly measures the JT distortion; the splitting between the (311)/(131) and (202)/(022) pairs of 
Bragg reflections is proportional to the order parameter. Our experimental results are compared 
to mean field calculations, taking into account all possible orientations of the grains relative to the 
applied field, and qualitative agreement is obtained. 

PACS numbers: 75.25.+z,75.10.-b 
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Phase transitions induced by quadrupolar interactions have recently received much at- 
tention, mostly due to their influence on the magnetic properties of the sample, e.g. in the 
giant-magneto-resistive effect of manganite compounds 



multi-ferroic materials 



l| and the magneto-electric effect in 



l3|. The cooperative Jahn- Teller (JT) effect arises from the same 
quadrupolar interactions 4j, and may indeed be interpreted as ferro-quadrupolar order. In 
JT compounds, the balance between magnetic and quadrupolar effects can be tuned by 
varying the strength of an externally applied magnetic field. To date, however, only a small 
amount of experimental evidence has been reported Q| , as the required magnetic fields tend 
to be rather large. 

In this letter we present the first direct experimental evidence for the modification of the 

both qualitative 



JT effect by a external magnetic field. This effect was predicted by 
and quantitative theories, but previous experiments 
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12] observed only 



indirect signatures of this behavior. 

Terbium ortho- vanadate, TbVC>4, along with DyVC>4, is a textbook example for a JT 
transition Q] induced by quadrupolar interactions between the Tb 4/ moments, mediated 
through phonons. At high temperatures, TbVC>4 crystallizes in the tetragonal zircon struc- 
ture with space group IAi/amd with lattice parameters a t — b t — 7.1831(3) A and 
c = 6.3310(4) A. Upon lowering the temperature through Tq ps 33 K, it undergoes a co- 
operative JT transition: The crystal spontaneously distorts along the [110] direction (B 2g 
strain, ^-symmetry distortion) to the orthorhombic space group Fddd with lattice param- 
eters a Q = 10.239(2) A, b Q = 10.029(2) A, and c = 6.3154(13) A at 22 K. The distortion 
is surprisingly large, reaching e 5 = 2(a a — b )/(a + b a ) = 2% at 22 K [l^ . and increasing 
further toward lower temperatures. 

The JT transition has been studied extensively, both experimentally and theoretically. 
The early studies, however, focused on the sample's properties in the absence of applied 



magnetic fields 4J. The effect of large externa 
recently. Quantitative mean-field calculations 



^magnetic fields on TbVC>4 was studied only 
found that the JT distortion is suppressed 



when fields above ~ 29 T are applied along the c-axis of the sample, in good agreement with 
susceptibility measurements 9J. Diffraction experiments have so far not been carried out, as 
the required combination of high magnetic fields and X-ray diffraction equipment has only 
recently become available jl4 . 115I. [if]] . 

The experiments were performed on the DUBBLE CRG beamline, BM26B, at the ESRF 
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TABLE I: Bragg reflections of TDVO4, scattering angle at E = 21 keV, and the splitting due to 



the JT-effect. Values are based in the lattice parameters given in Ref. Il3| for T = 22 K. 
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(HKL) 


29 


5(26) /e s 




[deg] 


[deg/%] 


(111) 


7.14 





(220) 


9.45 





(311)/(131) 


11.78/11.94 


0.076 


(202)/(022) 


12.61/12.69 


0.036 


(400)/(040) 


13.24/13.52 


0.134 


(222) 


14.32 





(331) 


15.19 





(113) 


16.81 






171 ]. using; the experimental configuration described in [16| . Flux-grown single crystals of 



TWO, 



18 



Q. 



19] were ground into a fine powder and embedded in low molecular weight 



polyvinylpyrolidone in order to prevent movement of the powder grains due to magnetic 
forces, while at the same time improving the thermal contact. The sample was then mounted 



in a pulsed magnet and cryostat assembly 



16J. 



The data shown here were taken by accumulating rs 45 magnetic field pulses per spectrum. 
For each field pulse a mechanical shutter exposed the image plate detector for 4.9 ms centered 
around the maximum field. Representative powder diffraction spectra are shown in Fig. [TJ 

The JT transition manifests itself as splitting of some, but not all, powder lines (see 
tabled]). For small distortions, C 1, the orthorhombic lattice parameters can be approx- 



imated as a Q ~ (1 + he 5 )a and b Q ~ (1 — he s )a, where a = h(a Q + b Q ) « \^2a t . To first order 



in e 5 , the splitting of a pair of Bragg reflections (H, K, L) Q and (K, H, L) Q is then given by 



5(26) 



A 2 K< 



H 2 



a 2 sin(2#) 



tan(0) 



K 2 -H 2 



H 2 + K 2 + ^L 2 



(1) 



Our experimental configuration, with photon energy E — 21 keV (A = 0.59 A), allowed us 
to observe scattering angles up to 29 = 12.8°, covering the reflections, in the orthorhombic 
notation, (111), (220), (311), (131), (202), and (022). Within this range, only the (311)/(131) 
and the (202) / (022) pairs are sensitive to the JT distortion. 
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FIG. 1: X-ray powder diffraction spectra of TbV04. The inset shows the (311)/(131) and 
(202)/(022) pairs of reflections that are sensitive to the orthorhombic distortion, e s . 

In samples with large magneto-crystalline (MC) anisotropy, such as TbV04, the phase 
diagram depends not only on the magnitude of an applied magnetic field, but also on its 
orientation relative to the crystal symmetry axes. The grains in a powder sample are aligned 
randomly, so that for a macroscopic measurement, e.g. magnetic susceptibility, the direction 
of the field has to be averaged over the full solid angle of 4ir. 

In a diffraction experiment, however, the orientation is constrained: Any grain contribut- 
ing to the diffraction signal at the angle 29 must be oriented such that the Bragg planes 
of the corresponding reflection (H, K, L) form the angle 9 with the incident beam. The 
powder averaging must then be taken only over the azimuthal angle, i/j, i.e. rotation about 
the scattering vector Q = (H, K, L). 

In our experiment the external field, [1qH, was aligned along the incident beam axis; the 
angle between the field and the Bragg planes is therefore the Bragg angle, 9. 

For systems with substantial MC anisotropy any simulation of the powder diffraction 
spectrum must thus be carried out peak by peak: For each Bragg reflection the magnetic field 
direction and the corresponding influence on the phase diagram and the lattice parameters 
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must be calculated as function of if>, yielding, in this case the distortion e s (if>), and the 
corresponding scattering angle, 29 ca \ c (HKL,e s (ip)). The resulting diffraction patterns are 
then averaged over the azimuthal angle, taking into account the instrumental resolution. 

In order to quantitatively describe the effect outlined above for TbV04 we have performed 
comprehensive mean field calculations of the magnetoelastic distortion of TbV04 as function 
of the strength and direction of the externally applied magnetic field. Our calculations closely 
follow those of j§| and and will be published in full detail elsewhere. In particular, we 
have used crystal field parameter set 2 of Q], and a quadrupolar constant G s = 130 mK. 

The order parameter (OP), e 5 , is proportional to the expectation value of the quadrupole 



operator 



22], Qxy {Pxy) j 



B 
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6 ~ qsQw (2) 

where P xy = \{J x J y + J y J x ). The coefficient B s /C% = 13.9 • 1(T 4 was determined from the 
zero-field, low temperature value of the spontaneous deformation, e s . 

We first consider the influence of the applied magnetic field on the quadrupole moment 
Q xy . Two effects must be taken into account: Changes in the magnitude of the distortion 
will result in variations of the splitting between the pairs of peaks, whereas (for a given 
grain) the preferential population of domains with positive or negative distortion will lead 
to a magnetically induced texture that manifests itself in a shift of spectral weight from one 
partner to the other. 

A magnetic field along the c axis leads to a non-vanishing expectation value of (J z ) in 
competition with the OP. A magnetic field along a Q or b a , on the other hand, will induce a 
magnetization (J x ± J y ) along the corresponding axis. Above the JT transition a magne- 
tostrictive orthorhombic distortion of 5-symmetry is induced. In the ordered state, T < Tq, 
however, one has to distinguish between a field applied along a Q and b Q \ Along the two axes, 
the susceptibility will be of slightly different value, thus breaking the degeneracy of ±e s . 
Crossing the phase boundary in the presence of such a magnetic field will then lead to the 
preferential population of one of these domains. Fields along the [110] o directions should 
lead to a distortion with 7-symmetry (Bi 9 -type), as observed in DyV04 20|, |2l|]. Finally, 
if the external field is applied along an arbitrary direction, then both S and 7 strains are 
induced simultaneously, lowering the symmetry of the crystal to monoclinic. Fortunately, 
the 7 magnetostrictive distortion in TbV04 is two orders of magnitude smaller than the S 
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FIG. 2: Calculated dependence of the quadrupole moment Q xy on the external magnetic field for 
different azimuthal angles in the interval (0-100°) at the diffraction peak (311). Left: T = 7.5 K, 
Right: T = 39 K. 

distortion and can therefore be neglected. 

Typical theoretical results for the (311) diffraction peak at T = 7.5 K and 39 K are shown 
in Fig. [2j The (311) partner of the (311)/ (131) pair corresponds to a positive quadrupole 
moment, Q xy > (see eq. [T|). 

At T = 39 K > Tq there is no spontaneous quadrupole moment at H = 0. For azimuthal 
angles ip below ~ 22° a positive quadrupole moment is induced, whereas for larger azimuthal 
angles the induced quadrupole moment is negative, i.e. corresponding to the (131) partner. 
We thus expect the magnetic field to induce a magnetoelastic splitting with unequal domain 
populations of ~ 33% Q xy > (contributing to the (311) diffraction peak) and ~ 67% 
Qxy < (contributing to the (131) diffraction peak) when averaging over the full interval 
(0-360°). 

At T = 7.5 K < Tq the situation is different: A sizable spontaneous quadrupole moment 
is present, and in a zero-field cooled sample domains with positive and negative Q xy will 
be populated almost equally. Upon applying an external field the quadrupole moment is 
reduced, and the degeneracy of ±Q xy is lifted, so that (depending on their azimuthal angle) 
domains may invert their quadrupole moment. This transition is first order and therefore 
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FIG. 3: Comparison between calculated and measured spectra for different temperatures and fields. 
Left: T = 7.5 K. Right: T = 39 K. At T = 39 K and 30 T, the dotted curve (in red) corresponds 
to the calculations including the magneto caloric effect, and the green line to the calculations with 
the coefficient B s /C$ reduced by 25%. 

subject to pinnin g on defects, grain boundaries, surfaces, and elastic interactions between 



twinned domains 



13|- 



The calculated dependences Q X y{H) for fixed T and different azimuthal angles ip allow 
us to simulate the X-ray powder diffraction data. For the (202) / (022) pair the values were 
averaged over the interval 0° < ip < 180° since the two-fold symmetry persists. For the 
(311) /(131) pair the averaging was performed over 0° < ip < 360°. The parameters of the 
resolution function were fitted to the experimental data at T > Tq and H = 0. Differ- 
ent parameter sets were used for the (311) / (131) and (202)/(022) pairs, and the intensity 
parameter had to be adjusted to match the data at low temperature and high field. 
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FIG. 4: Average order parameter as function of applied magnetic field for selected temperatures. 
Continuous lines: Experimental data. Dashed lines: Theory 

At T > Tq, calculations for the (202)/(022) pair adequately reproduce the almost sym- 
metric splitting of the line by the magnetic field observed in the data (Fig. [3]). The calculated 
magnitude of the splitting, however, is ~ 25% larger than the experimental one (Fig. [4]). 
The origin of this discrepancy remains unclear. In particular, it does not seem to be related 
to the value of the interaction parameters used in the calculations since the value of the 
spontaneous splitting at H = is correctly described (Fig. H]). Similar deviations between 
theory and experiment are observed at the (311) / (131) pair. Calculations describe the shape 
of the (311) line in different fields: the intensity ratio for two components of the split peak 
is about 1:2. But again the calculated value of the splitting is larger than the experimental 
one by about 25% (see Fig. [3]). The OP may be decreased by sample heating, e.g. due to 
the magneto-caloric effect. A reasonable temperature increase of about 8-10 K, however, 
cannot fully explain the observed discrepancy. 

The experimental observations at T > Tq and H = 30 T are well described by our 
calculations when the coefficient B 5 /Cq is reduced by 25%; this value, however, is too small 
for the spontaneous deformation. This may be interpreted as the variation of the quadrupolar 
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constant G 5 under high magnetic field, as was observed for the quadrupolar constant G 7 of 



22 1. However 



DyVC>4 and considered within an improved "compressible Ising model'' 
for TbV04 no dependence of G s was observed when the OP increases as the temperature 
decreases, and its low field behavior is described adequately by the mean field theory 4j. 

For the quadrupole ordered phase the results of the numerical calculations for the (202) 
and (311) peaks are compared with the experimental data in Fig. [3j For both peaks the 
averaging were performed for two domains (e s and — e s ) for < tp < 360°, in steps of 
Ai[) = 5°. For the (202)/ (022) pair the calculations do not reveal a noticeable decrease 
of the splitting under external field of 30 T, while for the (311) / (131) pair a small shift is 
expected for the low-angle component. 

We examined, and refuted, several possible explanations for the discrepancy between 
the theoretical and experimental results. Firstly we ascertained that various weak pair 
interactions (bilinear, quadrupole of 7-, a- and e- symmetry) which were omitted in our 
calculations do not change the results noticeably. Then we estimated the effect of the sample 
heating due to the magnetocaloric effect in our experiment. In the adiabatic regime, the 
temperature change is large enough (about 25 K at 7.5 K) for the magnetic field orientation 
close to the easy magnetization axis. Additional studies are necessary to elucidate the 
influence of the magnetocaloric effect in our experiments. 

In summary, we have directly observed for the first time the effect of magnetic fields 
on the Jahn- Teller distortion of TbVC>4, and compared the experimental data to a detailed 
mean field theory. Due to the polycrystalline nature of the sample and its strong MC 
anisotropy the calculated spectra had to be averaged over all possible orientations of the 
powder grains relative to the applied magnetic field. The applied magnetic field was found 
to influence both the magnitude of the OP, as observed in the splitting of the (311) / (131) 
and (202)/(022) pairs of Bragg peaks, and the relative domain populations, reflected in the 
intensity ratio between the partners of a pair. Our theory is in qualitative agreement with 
the experimental results, even though small quantitative discrepancies persist. 

This new technique was found to be a powerful tool to study JT systems, and allows us 
to observe new aspects of magnetic behavior of the classical and well studied JT compound 
TbV04. These high magnetic field data may be used both to revise known theoretical 
models and to develop new, improved ones. 
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